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Abstract 

In this paper the nonlinear matrix equation X - A*X~ P A = Q with p > is investigated. We 
consider two cases of this equation: the case p > 1 and the case < p < 1. In the case 
p > 1, a new sufficient condition for the existence of a unique positive definite solution for the 
matrix equation is obtained. A perturbation estimate for the positive definite solution is derived. 
Explicit expressions of the condition number for the positive definite solution are given. In the 
case < p < 1, a new sharper perturbation bound for the unique positive definite solution 
is evaluated. A new backward error of an approximate solution to the unique positive definite 
solution is obtained. The theoretical results are illustrated by numerical examples. 
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1. Introduction 

In this paper we consider the Hermitian positive definite solution of the nonlinear matrix 
equation 

X-A*X- p A = Q, (1.1) 

where A, Q and X are n x n complex matrices, Q is a positive definite matrix and p > 0. This 
type of nonlinear matrix equations arises in the analysis of ladder networks, the d yna mic pro - 



gramming, control theory, stochastic filtering, statistics and many applications [ 1 

In the last few years, Eq. dl.U was inv estigate d in some special cases. For the nonlinear matrix 
ecmations X - A*X~ X A = O |ilj fljjjll fl9[l23|l. X - A*X~ 2 A = Q JHH, X - A*X"A = Q 
|E [3 and X s - A*X-'A = QU^ there were many contributions in the literature to the 
solvability, numeri cal solutions and p erturbation analysis. In addition, the similar equations 

MUllS, X+A*X- 2 A = Q JMIIllii, X+A*X-"A = Q 

Q 



x+a*x~ i a = q JlIMIlElGlEi 

I20I1 were studied by many scholars 



M IHl], X+A*X-iA 



QtM\ZA\m,X+A*X"A 
Q JTHSEl and X+ £ A*X- l Ai 
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In Q, a sufficient condition for the equation X - A*X P A — Q (0 < p < 1) to have a 
unique positive definite solution was provided. When the coefficient matrix A is nonsingular, 
several sufficient conditions for the equation X - A*X~ q A — Q (q > 1) to have a unique positive 
definite solution were given in B33I1 . When the coefficient matrix A is an arbitrary complex matrix, 
necessary conditions and sufficient conditions for the existence of positive definite solutions for 



the equation X-A*X~ q A — Q (q > 1 ) were derived in [ 3 1 ] . Li and Zhang in 12411 proved that there 
always exists a unique positive definite solution to the equation X - A*X~ P A = Q (0 < p < 1). 
They also obtained a perturbation bound and a backward error of an approximate solution for the 
unique solution of the equation X - A*X~ P A — Q (0 < p < 1). 

As a continuation of the previous results, the rest of the paper is organized as follows. Section 
2 gives some preliminary lemmas that will be needed to develop this work. In Section 3, a new 
sufficient condition for Eq. dl.ll ) with p > 1 existing a unique positive definite solution is derived. 
In Section 4, a perturbation bound for the positive definite solution to Eq. dl.ll ) with p > 1 is 
given. In Section 5, applying the integral representation of matrix function, we also discuss 
the explicit expressions of condition number for the positive definite solution to Eq. dl.ll ) with 
p > 1 . Furthermore, in Section 6, a new sharper perturbation bound for the unique positive 
definite solution to Eq. dl.lt with < p < 1 is evaluated. In Section 7, a new backward error 
of an approximate solution to Eq. dl.ll ) with < p < 1 is obtained. Finally, several numerical 
examples are presented in Section 8. 

We denote by C" x " the set ofnxn complex matrices, by f{ nxn the set of n x n Hermitian 
matrices, by / the identity matrix, by || • || the spectral norm, by || ■ \\f the Frobenius norm and 
by A m . ix (M) and A m \ n (M) the maximal and minimal eigenvalues of M, respectively. For A = 
(ai,..., a„) = (a/j) e C" x " and a matrix B,A®B = is a Kronecker product, and vecA is a 

vector defined by vecA = (a[, a T n ) T . For X, Y e <H nxn , we write X > F(resp. X > Y) if X - Y 
is Hermitian positive semi-definite (resp. definite). Let k = /i max (A*A), k = A m i n (A*A). 

2. Priliminaries 

In this section, we will give some preliminary lemmas that will be needed to develop this 
work. 



Lemma 2.1. A241 For every positive definite matrix X e <H nxn i ifO < p < 1, then 

x - P = ^n AI+X) - lA - PdZ 

n Jo 



(if) X-p = ^IH. I (AI + Xy l X(AI + Xy l A- p dA. 

PK J() 



Lemma 2.2. 4241/ There exists a unique positive definite solution X of X — A*X P A — Q (0 < 
p < 1) and the iteration 

X >0, X n = Q + A'X'^A, n = l,2,--- (2.1) 

converges to X. 

Lemma 2.3. [Ml 

(f) IfX e ( H" X ", then \\ e - x \\ = e-^W 
(ff) IfX € < H nXn and r > 0, then X~ r = ^ J™ e^s^ds. 

(ii) If A, B € C" x ", Then e A+B - e A = jT 1 e^ k Be t( - k+B) dt . 
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3. A sufficient condition for the existence of a unique solution of X — A*X p A=Q{p>\) 

In this section, we derive a new sufficient condition for the existence of a unique solution of 
X - A*X~ P A - Q (p > 1) beginning with the lemma. 

Lemma 3.1. &l\llf 

p>(plc)^, (3.1) 

then Ea. (ll.l\) has a unique positive definite solution X e \J3I, al], where a and /3 are respectively 
positive solutions of the following equations 



(X - A max (Q)) [XminiQ) + —) = K 



and 



(X - A min (Q)) | AmaxiQ) + — I =K- 



Furthermore, 
Theorem 3.2. // 

(( P k)^ -A min (Q)) 



A mi „(Q) </3<a. 



A m a.x(Q) + 



then Ea. il.lt has a unique positive definite solution. 
Proof. We first prove 

i 

ft > (pT) . 

Let 



^_ A max (Q)(A min (Q) P y 

< K < K < : , 

(P ~ D" +1 



fix) = (x- A min (Q)) | AmaxiQ) + — I - K. 



By computaiton, we obtain 



Define that 



xP 



f'(x) = —\A max (Q) + — 



- \p-i 

K 1 



X p 



Amax(Q) „ 2i -1,1 : 

x p + p A min (Q)x +l-p 

K 



g(x) = x p + p z A min (Q)x + 1 - p l . 

K 



(3.2) 



(3.3) 



Then g(x) is decreasing on [0, ( 'V"^* )^' ] ana " increasing on [( ^J"^^" )"*' , +°°), which im- 
plies that 



■,(Q)pA 7 



AmaAQ) J 



= (1 +P) 



(A min (Q)p)- Acm 



+ l- P 



According to the condition k < ^ss^L^^El t ^ follows that g m m > 0. Noting that 



which implies that f(x) is increasing on (0, +00). Considering the condition (13.3b . one sees that 

1 . 1 

f{{pi<) ''+' ) < 0. Combining that and the definition of /} in Lemma DTll we obtain /3 > (pF) . By 

Lemma [3~T1 Eq. dl.U has a unique positive definite solution. □ 



4. Perturbation bound for X - A*X~ P A = Q(p>\) 



Li and Zhang in 12411 proved that there always exists a unique positive definite solution to the 
equation X - A*X~ P A — Q (0 < p < 1). They also obtained a perturbation bound for the unique 
solution. But their approaches will become invalid for the case of p > I. Since the equation 
X - A*X~ P A = Q(p > 1) does not always have a unique positive definite solution, there are two 
difficulties for perturbation analysis to the equation X - A*X~ P A = Q (p > 1). One difficulty 
is how to find some reasonable restrictions on the coefficient matrices of perturbed equation 
ensuring this equation has a unique positive definite solution. The other difficulty is how to find 
an expression of AX which is easy to handle. 

Assume that the coefficient matrix A is perturbed to A = AA + A. Let X = AX + X with 
AX e <H nxn satisfying the perturbed equation 

X-A*X- p A = Q, p>\. (4.1) 

In the following, we derive a perturbation estimate for the positive definite solution to the 
matrix equation X - A*X~ P A — Q (p > 1) beginning with the lemma. 



Lemma 4.1. h31\] If 

p\\Af<A p :'{Q), 



then Eq. il.lt has a unique positive definite solution X, where X > A m i n (Q)I. 
Theorem 4.2. // 



mil < y^H^ and HAAII < ^| — _ m< (4 . 2 ) 



then 

X-A*X- p A = Q and X-A*X- p A = Q 
have unique positive definite solutions X and X, respectively. Furthermore, 

\\X-X\\ ^ {2\\A\\ + HAAII 



Pll A p :'(Q)-p\\Af 



\\AA\\=g. 



Proof. By d4"5t . it follows that ||A|| < ||A|| + ||AA|| < J^^. According to Lemma gj] the 

condition (I4.2l i ensures that Eq. dl.U and Eq.( l4.U have unique positive definite solutions X and 
X, respectively. Furthermore, we obtain that 

X>A mm (Q)I, X>A min (Q)I. (4.3) 

Subtracting (14. U from (11. Il l gives 

AX = A*X- p A- A*X- p A = A\X~ P - X~ P )A + AA*X~ p A + A*X _P AA. (4.4) 
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By Lemma l23l and inequalities in (14.3b . we have 



\\AX + A*X- p A~A*X- p A\\ 

i r°° 

= \\AX + A*— {e- sX -e- sX )s p - l dsA\\ 

np) Jo 

= \\AX + A*-^- [ [ e- (1 -' )sX AXe- ,sX dts p dsA\\ 

Jo Jo 

> ||AX|| - l|A||2||AX|1 f M f' \\e-^ x \\\\e-' sX \\dts p ds 

r(/>) Jo Jo 

> iiaxii - ||A||2||Ax|1 r r 1 e -(n S i-(i) e -^m tfrfs 

Hp) Jo Jo 

> HAXII - l|A||2|l y 11 r C e -(l-')sA min (Q) e -,sA min (Q) dts P ds 

Hp) Jo Jo 

= ||AX||- I|A||2||AZ " f~ fe-^dtsQds 
T{p) Jo Jo 



r(p + l) ||A|| 2 ||AX|| 



IIAX|| 



mm v / 

^ p+l (Q)- p\\A\\ 2 

— ' IIAXII- (4.5) 

mm 



Noting fl4.2b , we have 



^T(0-rf>o. 



Combining ( 14.41 ) and ( 14.5b . one sees that 

&Q)-P\\A\\- 



\\AX\\ < ||AA*X- p A+A*X- p AA||<(||AA||+2||A||)||AA||||X- p || 



A P+ \Q) 

mm x 

< (||AA||+2||A||)||AA||^(Q). 



which implies that 



IIAXH (||AA||+2||A||) 

in,,, ^ — — : - AA . 

M ^(0-pHAlP 

□ 



5. Condition number for X - A*X~ P A = Q (p > 1) 

A condition number is a measurement of the sensitivity of the positive definite stabilizing 
solutions to small changes in the coefficient matrices. In this section, we apply the theory of 
condition number developed by Rice 112811 to derive explicit expressions of the condition number 
for the matrix equation X - A*X~ P A = Q (p > 1). 

5 



Here we consider the perturbed equation 



X-A*X~ P A = Q, p>l, (5.1) 

where A and Q are small perturbations of A and Q in Eg. ( II. Il l, respectively. 

Suppose that p\\A\\ 2 < A P ^(Q) and p\\A\\ 2 < /^f(Q). According to Lemma|4J] Eq.(EB 
and Eq. (15.1l l have unique positive definite solutions X and X, respectively. Let AX — X - X, 
AQ-Q-Q and AA - A - A. 

Subtracting ( 15.11 ) from ( 11.11 ) gives 



AX = A*X" P A - A*X~ p A + AQ = A*(JT P - Z _P )A + AA*X _P A + A*X~ P AA + AQ 

I r°° ~ _ 

= -A* (e~ sX -e~ sX )s p - 1 dsA + AA*X~ P A+ A*X~ P AA + AQ 

np)Jo 

= -A* f f - X)e~' sX dts p dsA + AA*X~ P A + A*X~ P AA + AQ 

r(p) Jo Jo 

= -A*— ( e - (l -' )sX - e- {l - ,)sX )AXe-' sX dts p dsA + AQ 

r(p) Jo Jo 

-A*— f f e - (x -' )sX AXe-' sX dts p dsA- (A* X- p AA- A* (X + AXy p AA) 
_T(p)J J„ 

+A*X~ P AA - (AA*X- p A - AA*(X + AX)~ P A) + AA*X~ P A 

= A*— f f f e- {X - m){l -' )sX AXe- m{X - ,)sX AXe- tsX dm{\-t)dts p+l dsA + AQ 
Jo Jo Jo 

-A*— f f e~ ( ' 1 ~ i)sX hXe~ t! ' x dts l 'dsA + AA*X~ P AA +A*X~ P AA + AA*X~ P A 
Hp) Jo Jo 

-A* — e^-'W^AXe^W^AA 
Hp) Jo Jo 

r(p)J Jo 

Therefore 

AX+A*— f f e- (1 - ,)sZ AXe- teZ dw p rfiA = E + h(AX), (5.2) 
T(p) Jo Jo 

where 



B = X~ P A, 

E = AQ + (B*AA + AA*B) + AA*X~ P AA, 

h(AX)=A*— f f f e- (l - m)(l - ,)sX AXe- m{X -' )sX AXe-' sX dm{\-t)dts p+l dsA 
T(p) Jo Jo Jo 

-A*— r C e^- t)s(X+Ax ^AXe-' sX dts p dsAA 
T(p)Jo Jo 

-AA*— r f' e- (l ->^ x+AX >AXe-' sX dts p dsA. 
T(p)J Jo 

6 



Lemma 5.1. If 



P\\A\\ 2 < A^(Q), (5.3) 



then the linear operator V : 7Y /IXrt — » 7Y" XM defined by 

i poo 

VW = W+ A*e- a - t)sX We-' sX Adts p ds, We<H nxn . (5.4) 

Hp) Jo Jo 

is invertible. 

Proof. Define the operator R : <H nxn -> "H**" by 



RZ=— f f A , «- ,1 -" ,A 'Z(-" v .Wh%.Zr//" 

r(p)J Jo 



it follows that 

Then V is invertible if and only if I + R is invertible. 

According to Lemma [2731 and the condition (15.31 1. we have 



i r°° r 

HRWII < ||A|| 2 ||W|| — \\e-^ x \\\\e-' sX \\dt S p ds 

n>) Jo Jo 

= IIA|| 2 ||W||||-i- f (V A V - :V ',//vV, 

r(p) Jo Jo 

< l|A|| 2 ||W||||-i- r [ e^-^^e-^^dtsPds 
r(p)J Jo 



i r°° 

l|A|| 2 ||W||— e -^«V^ 
n» Jo 



mf \\w\\ < \\w\\, 



A p+1 (0) 



which implies that ||R|| < 1 and I + R is invertible. Therefore, the operator V is invertible. □ 
Thus, we can rewrite (15.2b as 

AX = V~'Ag + V~'(B*AA + AA*B) + Y~ l (AA*X~ p AA) + V _1 (/z(AX)). 

Obviously, 

AX = V^AQ + V- l (B*AA + AA*B) + 0(\\(AA, A0|| 2 ), (AA, AQ) 0. (5.5) 

By the theory of condition number developed by Rice [6], we define the condition number of 
the Hermitian positive definite solution X to the matrix equation X - A*X~ P A — Q (p > 1) by 

c(X) = lim sup , (5.6) 

S ^°\\(f,f)\\F<s $ d 

where 77 and p are positive parameters. Taking £ = 77 = p — 1 in ( 15.6b gives the absolute 
condition number c a b s (X), and taking £ = 77 = ||A||f and p = HQIIf in ( 15.6b gives the relative 
condition number c re i(X). 
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Substituting ( 15.5b into ( 15.61 ), we get 



c(X) = - 



max 

' AA A£K 



!|V~'(A(2 + ZTAA + AA*B)||f 



i 

- max 

E eC" x ",H e<H n * 



Y- l (pH + T](B*E + E*B))\\ F 
\\(E,H)\\f ' 



Let V be the matrix representation of the linear operator V. Then it is easy to see that 

V = I®I+— [ f (e- tsX Af ®(A*e- (1 - ,)sX )dts p ds. (5.7) 

r(/>) Jo Jo 

Let 



V- 1 = 5 + zZ, 

v~' (/»#*) = y-'(/®(x~W) = Ui + ia u 

V~ l (B T = V~\{X~ p Af ®/)n = E/ 2 + iflfe, 



(5.8) 



5 -2 
Z 5 



Qi + Q 2 U x - U 2 



(5.9) 



vec// = x + iy, vec/s = a + ifo, g - (x T ,y T ,a T ,b T ) T , 

where x,y,a,Z> e 7?" 2 , S\Z, C/i, I/ 2 ,Qi,Q 2 e 7?" 2x " 2 ,M = (E,H),i = V~f, n is the vec- 
permutation matrix, i.e., 



vec E T — II vec £. 



Furthermore, we obtain that 



1 WV-'ipH + rjiB'E + E'BMF 
c(X) = - max 

1 llpV-'veci/ + 7?y- I ((/ ® B*)\ecE + (B T ® 7)vec£*)|| 

= - max 

£ M*0 ||(vecE,vecff)|| 

1 ||p(5 +iZ)(x + iy) + 7 ? [(f/i +iQi)(fl + i^) + (J/ 2 +iQ 2 )(fl-i^)]|| 

I M^O ||(vec£,vec#)ll 

1 ||(p5 c ,^ c )g|| 
= - max 

? g #0 Ml 
= ~ ||(pS c , 1/^)11, EeC* n ,Hz<H nxn . 



Then we have the following theorem. 
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Theorem 5.2. // p||A|| 2 < (Q), then the condition number c(X) defined by ( 15. 61 ) has the 
explicit expression 

c(X) = I || (pS e , nU c )\\, (5.10) 



where the matrices S c and U c are defined by ( l5.7D -( T5T9l ). 
Remark 1. From (15.10b we have the relative condition number 



™ II (UQUfSc, \\a\\ f u c ) II ..... 

Crd(x) = m • (5 ' n) 



5.7. The real case 



In this subsection we consider the real case, i.e., all the coefficient matrices A, Q of the matrix 
equation X - A*X~ P A — Q (p > 1) are real. In such a case the corresponding solution X is also 
real. Completely similar arguments as in Theorem l5.2| give the following theorem. 



Theorem 5.3. Let A, Q be real, c(X) be the condition number defined by ( 15.61 ). If p||A|| < 
A^(Q), then c(X) has the explicit expression 

c(X) = | || (pS r , nU r )\\, 



where 



S r = (/®/+J_ f P \ e -' sX A) T ®{A*e- il -' )sX )dts<>ds\ 

\ n» Jo J,, / 



u r = s r \j® (A T x- p ) + ((A T X- p ) ® /)n]. 

Remark 2. In the real case the relative condition number is given by 

\\ F S r , \\A\\ F U r )\\ 



CreliX) = 



6. New perturbation bound for X - A*X~ P A = Q (0 < p < 1) 

Here we consider the perturbed equation 

X-A*X- p A = Q, 0<p<l, (6.1) 

where A and Q are small perturbations of A and Q in Eq. dl.U . respectively. We assume that X 
and X are the solutions of Eq. dl.U and Eq.( 16.11 l. respectively. Let AX — X - X, AQ — Q — Q and 
AA =A-A. 

In this section we develop a new perturbation bound for the solution of Eq. dl.lb which is 
sharper than that in Theorem 3.1 1I24I1 . 

Subtracting Eq. dl.U from Eq. d6.ll ). using Lemma lXTl we have 

AX + I [(AI + Xy l A]*AX[(AI + Xy l A]A- p dA = E + h(AX), (6.2) 

n Jo 

9 



where 



B = X~ p A, 

E = AQ + (B*AA + AA*B) + AA*X- p AA, (6.3) 

h(AX) = !!^A*f (AI+Xy l AX(AI+X+AX)- l AX(AI+X)- l A- p dAA 
n Jo 

- !!^A*f (AI+Xy l AX(AI+X+AX)- l A- p dAAA 
n Jo 



sin pn r 
— AA ( 

n Jo 



o 

(^/+X)" 1 AX(/l/+X+AX)" 1 ^ /, d/M. 



70 

By Lemma 5.1 in Q, the linear operator L : -K" x " -> "W"*" defined by 

LW = W + | [(/I/ + xy l A]*W[(AI + xy l A]A- p dA, W e <H n 

n Jo 

is invertible. 

We also define operator P : C" x " -> < W nx " by 

PZ = L -1 (.B*Z + Z*B), ZeC" x ", i = 1,2,- •• ,m. 



Thus,we can rewrite (16.21 > as 

AX = L~'Ag + PAA + L-\AA*X- p AA) + \T x Qi(AX)). (6.4) 

Define 

||L _1 ||= max IIL-'WH, ||P|| = max ||PZ||. 

W € ( H nxn Z € C" x " 

\\w\\ = l ||Z|| = l 

Now we denote 

/ = HL-'ir 1 , <r = irt ? = iix-'h, b = iipii, 77 = /^iiaii 2 

e = y||Aei| + «||AA|| + ^||AA|| 2 , cr = |#(2||A|| + ||AA||)||AA||. 



Theorem 6.1. // 

1(1 -of 

a < 1 ana! e < , (6.5) 

((I + Ict + 2t] + 2 y/(lcr + 77X77 + /)) 

P-X||< s /i. 

1(1 +£e-cr)+ V/ 2 (l +^e-o-) 2 -4/^e(/ + 77) 

Proo/ Let 

/(AX) = L^Ag + PAA +L"'(AA*^''AA) + L" 1 (MAX)). 

Obviously, / : "TY"*" — > <7Y" X " is continuous. The condition ( 16.51 ) ensures that the quadratic 
equation ((I + rj)x 2 - 1(1 + — cr)x + le — in x has two positive real roots. The smaller one is 

lie 

iU = 



1(1 + £e - cr) + ^P(l + le - o-) 2 - 41^(1 + 77) 
10 



Define Q. = {AX e <H nxn : \\AX\\ < p,}. Then for any AX e Q, by tf63T >. we have 

2e 



\\X- x AX\\<\\X-'\\\\AX\\<t^<t- 



\ + e -cr 



= 1 | ftr + o--l | -2(1 -<r)(Zo- + 17) ;1 

1+fe-cr" (7cr + Z + 2?7Xl +^e-cr) 

It follows that / - X~ l AX is nonsingular and 

\\I-X- l AX\\< ' ' 



Therefore 



1 £ , p , ||AX|| 2 

||/(AX)|| < -||Aei| + n||AA|| + ^||AA,|| 2 + ^||A|| 2 



/" ~ / i" 1 -f||AX|| 
P r ^(2||A|| + ||AA||)||AA|| ■ ilAA!! 



r 3V i-fllAX|| 

gjlAXjj 7?^||AX|| 2 

< € + Z 777— ~~~7~ + 



i-{\\ax\\ i(i -am) 

< EH 1 = 11, , 

l-f/i. «Wm.) 

for AX e Q. That is £ O. According to Schauder fixed point theorem, there exists AX, e D. 
such that /(AX,) = AX,. It follows that X + AX» is a Hermitian solution of Eq.( l6. lb . By Lemma 
2721 we know that the solution of Eq.d6TTl) is unique. Then AX, = X - X and \\X - X\\ < p„. □ 



7. New backward error for X - A*X~ P A = Q (0 < p < 1) 

In this section we evaluate a new backward error of an approximate solution to the unique 
solution, which is sharper than that in Theorem 4.1 112411 . 

Theorem 7.1. Let X > be an approximation to the solution X of di.il ). i/'||X _ iA|| 2 ||X _I || < 1 
and the residual R(X) = Q + A*X~ P A — X satisfies 



ITOII < min 1, l -=r , wtere ft = (1 - ||^- 2 A^HXT- 1 |IMmi„(X) + \\R{X)\\ > 0, (7.1) 

2 I 2i min (X)J 

?/zen _ 

\\X - X\\ < 6\\R(X)\\, where 6 = mM = (7.2) 

ft + ^6] - 4A~(X)\\R(X)\\ 

To prove the above theorem, we first verify the following lemma. 

Lemma 7.2. For every positive definite matrix X e •H" x ", < p < 1, ifX + AX > (l/v)7 > 0, 
then 

||A*((X + AX)-" - X-')A\\ < p (IIAXII + v|IAX|| 2 )||X-f AlPlIX" 1 !!. (7.3) 
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Proof. \fX + AX > (1/v)/ > 0, then 

\\A*((X + Axy p -X~ P )A\\ 



A * j~( (AI+x+ ax)-! _ {AI+ xyl j rP dA y 

Ua" (AI + xy l AX(A I + X)- l A- p dA A||j 
^ |||A* (AI + Xy l AX(AI + X + AXy l AX(A I + Xy l A- p dA A||j 



sin 

+ 



< p IIA^AIIIIX^HIIAZII + p IIA'X-MHvllAXH 2 !^- 1 ! 



p(\\AX\\ + v\\AX\\ 2 W^A\\ 2 \\X- l \ 



□ 



Proof. Let 

*P = {AX e Tf"*" : ||AX|| < e\\R(X)\\}. 
Obviously, ¥ is a nonempty bounded convex closed set. Let 

g(AX) = A*((X + AXy p - X- p )A + R(X). 

Evidently g : T h </y mx " is continuous. We will prove that §( V F) £ "F. For every AX e x ¥, we 
have 

AX > -6\\R(X)\\I. 

Hence 

X + AX >X- 0\\R(X)\\I > (A min (X) - e\\R(X)\\)I. 
Using ( 17. U and ( 17.2b . one sees that 

fllKftll = ^■WI^WH < 2^pgH < ^ 



* + 



- 4A min (X)\\R(X)\\ 



Therefore, (A^X) - 8\\R(X)\\)I > 0. 
According to ( 17.3l l. we obtain 

\\g(AX)\\ 



< p (HAZII + ^)l|3r5 AlhlX- 1 1| + \\R(X)\\ 

Amin(X) - WWW 



I t n\\T>lv\\\\l \ 



(e\\R(X)\\) 2 



6\\R(X)\\ + 

A min (X) - e\\R(X)\\ 

0\\R(X)\\. 



(p\\X-2A\\ 2 \\X- l \\) + \\R(X)\\ 



By Brouwer's fixed point theorem, there exists a AX e *P such that g(AX) - AX. Hence X + AX 
is a solution of Eq. fll.ll i. Moreover, by Lemma [2721 we know that the solution X of Eq. fll.U is 
unique. Then 

||X-X|| = ||AX||<0||/?(X)||. 

□ 
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8. Numerical Examples 

To illustrate the theoretical results of the previous sections, in this section four simple exam- 
ples are given, which were carried out using MATLAB 7.1. For the stopping criterion we take 
£k+i(X) = \\X k - A*X k P A - Q\\ < l.Oe - 10. 

Example 8.1. We consider the matrix equation 

X-A*X- h A = I, 



where 



A = 



\\A 



An = 



2 0.95 
1 



Suppose that the coefficient matrix A is perturbed to A — A + AA, where 

IQ-j 



AA = 



\\C T + C\\ 



(C T + C) 



and C is a random matrix generated by MATLAB function randn. 

We compare our own result & in Theorem \6.1w ith the perturbation bound f» proposed in 
Theorem 3.1 h24\l . 

The condition in Theorem 3.1 [24] is 



con\ = ^||A|| 2 +<r-||A||-||AA|| > 0. 
The conditions in Theorem \6.1\ are 

1(1 -o-) 2 



conl = 1 — o" > 0, con3 = 



£(/ + o-l + 2r] + 2 V(/cr + T])(i] + I)) 



-£>0. 



By computation, we list them in Table\l\ 



Table 1: Conditions for Example l8.1l with different values of j 



i 


4 


5 


6 


7 


conl 


0.0455 


0.0456 


0.0456 


0.0456 


con2 


0.9999 


1.0000 


1.0000 


1.0000 


con3 


0.3957 


0.3959 


0.3959 


0.3959 



The results listed in Table\l\show that the conditions in Theorem 3.1 h241 and Theorem \6.1\ 
are satisfied. 



By Theorem 3.1 in [24] and Theorem \6.1\ we can compute the relative perturbation bounds 
€*' IH' res P ect ' ve fy- These results averaged as the geometric mean of 10 randomly perturbed 
runs. Some results are listed in Table\2\ 

The results listed in Table\2\show that the perturbation bound given by Theorem \6.1\ is 



fairly sharp, while the bound ^ given by Theorem 3.1 in H24\] is conservati 



ve. 
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Table 2: Results for Example l8.1l with different values of j 



j 


4 


5 


6 


7 


ira 


6.8119 x 10~ 5 


4.2332 x 10~ 6 


4.3287 x 10~ 7 


5.5767 x 10~ 8 


f. 


2.6003 x 10~ 4 


2.1375 x 10~ 5 


1.9229 x 10" 6 


2.7300 x 10~ 7 


11X11 


8.8966 x 10~ 5 


6.5825 x 10- 6 


7.2867 x 10" 7 


9.3455 x 10~ 8 



Example 8.2. Consider the equation 

X - A*JT 3/4 A = Q, 

for 

0.2 -0.2 \ / 0.8939 0.2987 

0.1 0.1 j' Q ~ [ 0.1991 0.6614 

Choose Xo — 3Q. Let the approximate solution X k be given with the iterative method ( 12.71 ), where 
k is the iterative number. Assume that the solution X ofEq. fLl\ is unknown. 

We compare our own result with the backward error proposed in Theorem 4. 1 [24 ]. 

The residual R(X k ) = Q + A*X k P A — X k satisfies the conditions in Theorem 4.1 n24\l and in 
Theorem \7.1\ 

By Theorem 4.1 in[24] , we can compute the backward error bound 



X k - X \\< v,\\R(X k )\\, where v. = 



2\\X k \\\\X^\ 



i - 1 pr 5 Axr 1/2 ii 2 



4 IK*£ 

By Theorem YI .1\ we can compute the new backward error bound 

2/t m in(X:) 



||X* - X|| < 0||/?(X*)ll, where Q ■. 



6i = (i- pruiiwiWmtaWt) + WW 



h + - 4A min (X k )\\R(X k 



k "^k 

Let _ _ 

nWWM o\\R(x k ) 

K\ = — , K 2 = —=- 



\x k -x\\ ' \\x k -x\\ 

Some results are shown in Tabl$3\ 

From the results listed in Table\3\we see that the new backward error bound 6\\R(X k )\\ is 



sharper than the backward error bound v t \\R{X k )\\ in [24]. Moreover, we see that the backward 
error 0\\R(X)\\ for an approximate solution X seems to be independent of the conditioning of the 
solution X. 

Example 8.3. We consider the matrix equation 

X - A*X~ 3 A = 51, 
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Table 3: Results for Example l8.2l with different values of k 



7 



\\x k -x\\ 
y*\\R(x k )\\ 

e\\R(x k )\\ 

K2 



6.2131 x 10" 
2.5930 x 10" 
4.1734 

7.0053 x 10" 
1.1275 



1.5830 x 10" 
6.6257 x 10" 
4.1856 

1.7900 x 10" 
1.1308 



8.2486 x 10" 
3.5697 x 10" 
4.3277 

9.6440 x 10" 
1.1692 



6.0132 x 10" 10 
2.4646 x 10" 
4.0986 

6.6583 x 10"° 
1.1073 



where 

A , / 2 0.95 

A 



M = „ " T " „„ (C T + C) 



iiAoir \ o i 

We now consider the perturbation bounds for the solution X when the coefficient matrix A is 
perturbed to A — A + AA, where 

10" 

iic^+cir 

and C is a random matrix generated by MATLAB function randn. 
The conditions in Theorem \4.2\ are satisfied. 

By Theorem \4.2\ we can compute the relative perturbation bound g with different values of j. 
These results averaged as the geometric mean of 10 randomly perturbed runs. Some results are 
listed in Tabled 

Table 4: Results for Example 1 8 . 3 1 with different values of j 



j 


4 


5 


6 


7 


lix-xil 
\\x\\ 


1.1892 x 10" 


2.1101 x 10" 


2.4085 x 10" 


1.6847 x 10"° 


Q 


2.0791 x 10" 


3.5353 x 10" 


3.9573 x 10" 


3.2580 x 10"° 



The results listed in Table^show that the perturbation bound q given by Theorem \4.2\ is fairly 
sharp. 

Example 8.4. Consider the matrix equation X — A*X~ 3 A — Q, where 
. I 0.5 0.55 - 10" \ _ / 5 1 

A =l i i f e ii 5 

By Remark^ we can compute the relative condition number c re i(X). Some results are listed in 
Table\5\ 

The numerical results listed in the second line show that the unique positive definite solution 
X is well-conditioned. 
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Table 5: Results for Example l8.4l with different values of k 



k 


1 


3 


5 


7 


9 


Crel(X) 


1.2510 


1.0991 


1.0009 


1.0009 


1.0009 
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